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Abstract 

In six chosen scenarios for the dependence of the form factors involved in 
D'^ — > (j) transition, we have determined the allowed domain of x = A2(0)/Ai(0) 
and y = V{Q)/Ai{Q) from the experimentally measured ratios Rsi = ^{Df — > 
(j)l+ui)/T{D+ (/.7r+) and Rh = T{D+ (l)p+)/T{D+ (/)7r+) in a scheme that 
uses the Nc = 3 value of the phenomenological parameter ai and includes nonfac- 
torized contribution. We find that the experimentally measured values of x and y 
from semileptonic decays of favor solutions which have significant nonfactorized 
contribution, and, in particular, Rgi favors solutions in scenarios where Ai{q'^) is 
either fiat or decreasing with q^. 

PACS index: 13.25Hw, 14.40.Nd 
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I. Introduction 

The phenomenological parameters ai and 02, appearing in the description of charm- 
ed and bottom meson hadronic decays, are related, for Nc colors, to the perturbatively 
calculable Wilson coefficients Ci and C2 by fl], , 

ai,2 = Ci,2 + (1) 

In the description of two-body hadronic decays of D mesons, phenomenology seemed 
to support the view that Nc 00 limit applied ||^. This implied that ai(mc) ~ Ci(mc) ~ 
1.26 and a2{mc) ~ C2{mc) ~ —0.51 Q]. However, the same idea {N^ 00 limit) carried 
over to hadronic B decays failed. Nc —>■ 00 limit would have implied ai{mb) ~ 1.12 and 
a2(mc) ~ —0.26 @], yet two-body B decay data leave no doubt that a2{mb) is positive 

We have recently proposed that in i? ^ ■ip{ip{2S))+K{K*) decays, one ought to use 
Nc = 3 values of ai and 02, ai(mf,) ^ 1.03 and a2{mb) ~ 0.09, and absorb the nonfactorized 
contributions in the definition of effective, and process-dependent, Oi and 02 for processes 
such as B ^ PP, VP (P and V represent 0~ and 1~ mesons, respectively) which involve 
only a single Lorentz scalar. The idea of effective ai and 02 is not new; it has been proposed 
and discussed by various authors in the past p, ^ |10|, |ll|. For processes of kind B VV, 
involving three independent Lorentz structures (as three partial waves, S, P, and D, come 
into play), it is not possible to define effective ai and 02 0. Pursuing this idea, we 
had shown that color-suppressed decays B — * ip{ip{2S)) + K{K*) and the longitudinal 
polarization fraction in B —>■ + K* could be understood in all commonly used models 
provided one included a small nonfactorized amplitude, the amount depending on the 
model. This ameliorates the conundrum posed in [|1^ where it was shown that none of 
the commonly used models explained the polarization data in B ^ ip + K* decay within 
the factorization assumption. 

In the present paper we have carried this idea into the realm of charmed meson decays; 
in particular, we investigate the decays (f)TT'^,(f)p~^ and (pl^i^i. Experimentally, only 

relative rates are measured (following the notation introduced in [0), 

_ V{Dt ^ 0P+) 
r(D+^07r+)' 
_ T{Dt ^ 0/+^z) 
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Though the problem has been discussed in |T3[, we revisit it with the purpose of 
investigating the role of nonfactorized contribution in a formalism that uses the A''^ = 3 
values of ai{mc) and 02 (^c). The reason for choosing the Cabibbo-favored decays of D'^ 
is that the hadronic final states involve only a single isospin and, consequently, one does 
not have to worry about the interference effects due to final state interactions. 

To the best of our knowledge nonfactorization contributions were first discussed in 



charm decays by Deshpande et al. in 



In Section II, we describe the method and the calculation of the three decay rates out 
of which Rh and Rgi are constructed. The results are presented and discussed in Section 
III. 

II. Method and Calculations 
A. Definitions 

The effective Hamiltonian for Cabibbo-favored hadronic charm decays is given by 

= ^KsKTd {Ciiud)isc) + C2iuc)isd)} , (3) 

where (ud) etc. represent color-singlet (V - A) Dirac currents and Ci and C2 are the 
Wilson coefficients for which we adopt the following values, 

Ci = 1.26 ± 0.04, C2 = -0.51 ± 0.05. (4) 

The central values of Ci and C2 are taken from |Q , though the error assignments are ours. 
Fierz-transforming in color space with Ac = 3, we get 

1 1 ^ 

(uc) (sd) = - (ud) isc) + -J2 (^A"c) (sX^d) , (5) 

where A'^ are the Gell-Mann matrices. We define, for Ac = 3, 

C 

ai = Ci + — = 1.09 ±0.04, 
3 

C 

02 = C2 + = -0.09 ± 0.05. (6) 

o 



2 



We also define ^ the relevant hadronic matrix elements and form factors, 

(p+IMio) 

(vr+IMIO) 
msc)\Dt) 



e .q 



i\{mD^+m^)e*Ai{q ) - 



£ .q o 



2n , ^ 



+ 



(7) 



where = {pn^ — Pf/,)^- In addition, the following constraints apply at = 0, 

Ao{Q) = ^3(0), 

2m^A3{0) = {niD, + m^)Ai(O) - (m^^ - m^)A2{0). (8) 

With these preliminaries, we now study the rates for the processes fpl^vi-, (pTc^ and 



B. 



"z// decay 



In the notation of [|T3| 
dV 



X 



where, Hwit^) in each helicity state are defined by. 



1 + r 
2r 



with 



(l + r)2 



< < fl -r)^ 



m 



2 ' 



In the rest frame of Df, the magnitude of meson 3-momentum is given by 

Pit') 



with k{f) = [(l + r2-t2)2_4^2|l/2_ 



(9) 



(10) 



11^ 



121 
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C. — > 07r+ decay 



Using (3), (5) and (6), the decay amplitude for — > 07r+ is written as 

71 



A{Dt - 07r+) = ^V^sV:, {a, {<P'n+\{ud){sc)\Dt) + (07r+|i7(f)|Z}+)} , (13) 



where, from (5) 



H^^^ ^\Y.i^X''c){s\'^d). (14) 



We calculate the first term in the brackets of (13) in the factorization approxima- 
tion using the definitions in (6) and parametrize the second, the nonfactorized term for 
convenience of combining it with the first term, as follows, 

^^\H^:^\Dt) = 2m,Ue\pnAnml). (15) 



This nonfactorized parameter, ^4^"^, represents a spurion scattering H^^ + — > (f) + 7r~^, 
at the Mandelstam point s = rn^o^i t = and u = m^. 

Combining the factorized and the nonfactorized terms we obtain, 

A{Dj^(Pn+) = V2GpV,,V:,m^a'l^Uo{ml)e*.pD,, (16) 

where, 

. aM.^^}. (IT) 

In terms of the quantities defined in (10) - (12), the decay rate 07r~*" is given by 

in, 



327r \fnD 
where we have approximated / mj-,^ ~ and 



2\3 



fc(0)i/oo(0) = (l-0 



1 +r 
2r 



I + r 



2 



(19) 



4 



We choose not to work with the form (18) which involves Ao{0) in the definition of 
a^-^-^ but rather with a form in which Aq{0) is ehminated altogether in favor of ^i(O) and 
^42(0) via (8). In doing so, we obtain, 

r (Dj 07r+) = ^^\V,s\Wal (^]\{0)Hoo{0), (20) 



where. 



A;(0)^oo(0) = (1-r 



2\3 



1 + r 



2 



(21) 



2r 

D. Df (f)p'^ decay with a zero width p meson 

From (3), (5) and (6), the decay amplitude for Df — > is written as, 

A {Dt - = ^V^V:, {a, {cl>p+\{ud){sc)\Dt) + C, , (22) 

where H^^ is defined in (14). 

Again we calculate the first term in the brackets of (22) in the factorization approxima- 
tion and define the nonfactorized second term, for ease of combining it with the factorized 
first term, as follows, 

(0p+|H<f>|D,+) = im,f,[(mo.+mt)(e;.el)Ar(ml) 



2i 



-e,.,^e;%^f^y^V''^{ml)y (23) 



mo, + 

Further, because of lack of phase space in this decay mode, we retain nonfactorized con- 
tribution only to S- waves, that is, 

ArH) + 0, 

= ^^^K)=0. (24) 

The decay amplitudes in each helicity state are then given by 

G 

AmiP'^s = -^csyud^pfp{mD, + m^)ax {a^Ai(mJ) - M2(mJ)} , 

^±±W^#+) = ^K,T/>,/,(mz,.+m<^)ai{eAi(mJ)Tcy(mJ)}, (25) 
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where, with tp = mp/mr,^ 



m 



ml -ml 



2mpm^ 



2rtr. 



Ds 



k\tl) 



{mos + m^Ympm^ 2rtp{l + rf 



"^WA^Ds __ 

{mDs+m^y (l + r)2' 



(26) 



The decay rate can then be expressed in a form resembhng the expression given in 



[0 



r m: 



327r 



(27) 



with 



:28) 



E. (f)p^ decay with a finite width p meson. 

The finite width of p meson is taken into account by smearing the rates given in (27) 
and (28) over the p width by using a unit normahzed Breit-Wigner measure, p{t^)- This 



is accomphshed by the replacement [13, 15, 16] 

r-(l-r)" 



k{t')Hxx{e)p{t')dt' 



(29) 



with 



p{f)dt' 



m^ 



(30) 
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Two measures that have been used in the past are ^ 



TX 



it' - tiy + ijti 



(31) 



and 



No 



with 



4*2 



(32) 



(33) 



The appearance of A^^i 2 in (31) and (32) ensures that pi,2(^^) are unit normahzed. In (31), 
7p = rp/m^+, where Vp is the p width. In (32), 7(it^) is so chosen as to reflect the P-wave 
nature of the p meson 



lit' 



t2 



2 13/2 



tl - K 



(34) 



For a p width of 151.2 MeV and nip = 769.1 MeV, we found A^i = 1.0758 and N2 = 0.9946. 
In our calculation we have used the smearing function p2{t') with energy-dependent width. 

III. Results and Discussions 

In the results presented below we have used Vcs = Vud = 0.975 and only the central 
values of ai, a2 and C2 given in (4) and (6): ai = 1.09, a2 = -0.09 and C2 = -0.51. Defining 



X 



^2(0) 



y 



we can write the decay rate for </)7r+ as 

r(L)+ 07r+) = 0.2341 X lO^^l^i 

where, 



^ (1 - 0.3177X - 0.3192x^)^ s-\ 



(35) 



(36) 



We calculate T{D^ (pl^i^i) and the smeared T{Df in six different scenarios 

listed below. 

Scenario 1: All form factors extrapolate in as monopoles; Ailq"^) and ^2(0'^) with 
pole mass 2.53 GeV and V{q'^) with a pole mass 2.11 GeV (BSWI ). We get 

T{Dj ^ (pl+ui) = 0.1615 X 10^Vi(0)|^ (1 - 0.2840X + 0.0344x^ + 0.0144?/^) s~\ 
T{Dj 0P+) = 0.3972 X 10^Vi(0)p (^^ - 0.1493a;C + O.OllSa;^ + 0.0167^^) s-\ (38) 

Scenario 2: v4i(g^) extrapolates in as monopole with pole mass 2.53 GeV and ^2(9^) 
and V^q"^) as dipoles with pole masses 2.53 and 2.11 GeV, respectively (BSWII ). We 
obtain 

r(D+ ^ (pl+i^i) = 0.1615 X 10^^|Ai(0)p (1 - 0.2958X + 0.0367^^ + 0.0176?/^) 

T{D+ 0P+) = 0.3972 X 10^Vi(0)p (^^ - 0.1623x^ + 0.0137^^ + 0.0218^^) s'K (39) 

Scenario 3: Ai(q'^) flat, ^2(0'^) and V^q"^) extrapolate in as monopoles with pole 
masses 2.53 and 2.11 GeV, respectively. We obtain 

T{D+ (ptui) = 0.1423 X 10^2|Ai(0)p (l - 0.3098x + 0.0390^^ + 0.0163^^) 

r(D+ ^ 0P+) = 0.3297 X 10^Vi(0)|2 - 0.1656a;^ + 0.0142^^ + 0.02027/^) (40) 

Scenario 4: Ai(g^) decreasing linearly in g^, 

A(,^)=A.(0){l-(^)]. (41) 

and A2{q'^) and V(g^) extrapolate in g^ as monopoles with pole masses 2.53 and 2.11 GeV, 
respectively. This yields 

r(D+ ^ (pl+i^i) = 0.1262 X 10i2|^i(0)p (1 - 0.3361X + 0.0440^^ + OMSAy^) s-\ 
T{Dj 0P+) = 0.2742 X 10^Vi(0)|^ (^^ - 0.1835a;^ + 0.0171^^ + 0.0243?/^) s'K (42) 

Scenario 5: v4i(g^) and ^2(0'^) flat and V{q'^) extrapolates in g^ as monopole with pole 
mass 2.11 GeV Q. We get 

T{D+ (pl+ui) = 0.1423 X 10^Vi(0)p (l - 0.2980x + 0.0366x^ + 0.0163^^) 

T{D+ 0P+) = 0.3297 X 10^^\Ai{0)\^ (^^ - 0.1526a;C + 0.0122a;2 + 0.02027/^) (43) 
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(44) 



Scenario 6: Ai{q'^) decreasing linearly in according to [|13 

^.(.^)^^.(0){l-(^)], 

and v42(g^) and V^q"^) extrapolate in as monopoles with pole masses 2.53 and 2.11 GeV, 
respectively. This yields 

T{D+ (jyl^ui) = 0.1337 x 10^Vi(0)P (l - 0.3172x + 0.0403^^ + 0.0174?/^) s-\ 
T{Dj ^ 0P+) = 0.3000 X 10^Vi(0)|2 _ o.l667a;C + O.OMSx^ + 0.0222y2) s-\ (45) 

For data, we use 

Rf, = 1.86±0.26l°i9, [17] (46) 
Rsi = 0.54 ±0.10. [|18| (47) 

In the following, we discuss the analysis of the ratios Rgi and Rh separately. 

Rsi : We reiterate that we are using the Nc = 3 value of ai = 1.09 (central value only). 
Rsi is constructed from (36) for T{Df — > 07r+) and (38)-(45) for T{D^ 4>^'^^i) i^i various 
scenarios. The allowed region in (x, y) plane are shown in Fig.l for some selected scenarios 
as explained below. Firstly, we observe that no solutions were found for Xn = (see (37)) 
in scenarios 1 and 2 (BSWI and BSWII, respectively). In contrast, solutions were found 
in these scenarios in [jl3|; the difference lies in our use of Oi = 1.09 while Oi ~ 1.26 was 
used in |T^. Solutions, however, were found for = in all the other four scenarios. 
Secondly, the scenarios for Ai(g^) flat, and Ai(g^) and ^2(^7^) flat, for = or 7^ 
were almost indistinguishable, while those for BSWI and BSWII with Xn ^ were very 
similar. Consequently, we have chosen to plot only the results using scenarios 1, 3 and 4 
to keep Fig.l uncluttered. The experimental points are from Refs. [0, pO| and ||21[ and 
the plots are made for Xtt = and -0.5. 

The absolute rate, T{D'^ 07r+) and the branching ratio, B{D^ 07r+) can not be 
calculated in a model-independent way, but, for the record, in BSW model 0], one gets 

i?(D+^07r+) = 2.48% (X. = 0), 

= 4.03 % (x^ = -0.5). (48) 
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Particle Data Group list B{Df — > 07r+) = (3.5 ±0.4)%, though a direct measurement 
of the branching ratio does not exist. 

Clearly from Fig.l one notes that data prefer solutions with v4i(g^) flat or Ai(g^) 
decreasing with (scenarios 3 to 6). Further, solutions with a nonfactorized contribution, 
Xn = —0.5, fare much better than those with Xn = 0. In particular, E-687 data pO[ are 
consistent with all the six scenarios with Xn = —0.5. CLEO data [^| are consistent with 
scenarios 3 to 6 (Ai(g^) fiat or decreasing with g^) with x-k = —0.5. E-653 data |jT9[ do 
not admit a solution with —0.5 < Xtt ^ 0. 

Rh : The ratio Rh, eq.(2), is constructed from the rates T{D^ 07r+) given in (36) 
and T{Df 0P^)) for finite width p meson, given in (38) to (45). The allowed regions 
in (x, y) plane for the six scenarios are shown in Figs. 2, 3 and 4. We find that with 
Xtt = Xp = 0, solutions accommodate E-687 data pO[ in all scenarios, while CLEO data 
[ PD are accommodated in scenarios 4 and 6 (yli(g^) decreasing with g^). 

However with Xn = —0.5 and xp = 0.5, all three data points |T9|, ^ are accom- 
modated in all six scenarios. The allowed region in the {x, y) plane is now a band with 
an upper and lower branch. 

For the record, with Xp = 0.5, and using BSWII model [Q, one gets 

B{Df (f)p+) = 6.23 %. (49) 
Again, though this branching ratio is not directly measured. Particle Data Group |Q list 

B{Dt ct>P^) = (Q.btli) % (50) 

In conclusion, taking the Nc = 3 value of ai seriously, we have asked: What is the 
domain of x and y allowed by the ratios Rh and Rsi in six chosen scenarios for the 
dependence of the form factors, with and without nonfactorization contribution? And, 
what is the size and the sign of the nonfactorization contribution in order that the mea- 
sured values of x and y fall within the allowed domain obtained from Rh and Rsi'? We 
should emphasize that the experimental determination of x and y from semileptonic data 
is not model free as monopole extrapolation is assumed for all the form factors in data 
analysis. We find that an analysis with the inclusion of nonfactorized contribution fares 
much better in selecting an allowed domain of x and y consistent with the data. This 
is particularly true of Rsi where, in addition, the scenarios in which v4i(g^) is flat, or 
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decreasing with are favored over the scenarios where it rises with q^. As to the size 
and sign of the nonfactorized contribution, we have no explanation. 
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Figure Captions 



Fig.l: The regions enclosed by various curves and the axes are the domains of {x,y) 
allowed by Rsi (eqn. (47)) in different scenarios. The solid curve is for scenario 1 with 
Xtt = —0.5; the dashed (second from the right) and the dash-dashed (innermost) are 
for scenario 3 with = —0.5 and 0, respectively; the dotted (outermost) and the dot 
-dashed (second from left) are for scenario 4 with Xn = —0.5 and 0, respectively. The 
corresponding curves for scenario 6 (not shown to avoid cluttering) lie between the curves 
for scenarios 3 and 4. Also shown are the data points; A: E-687 |^D|, B: E-653 [|Ty] and 



C: CLEO PI 



Fig. 2: The domain (between two solid lines) of {x,y) allowed by Rh (eqn. (46)) with 
Xp = 0.5 and Xtt = —0.5. (a) in scenario 1; (b) in scenario 2. The region enclosed by 
the dashed lines and the axes both in (a) and (b) is the allowed domain of (x, y) with 
Xp = Xn = 0. Also shown are the data points; A: E-687 0, B: E-653 ^ and C: CLEO 

Fig. 3: Same as Fig. 2. (a) in scenario 3; (b) in scenario 5. 
Fig. 4: Same as Fig. 2. (a) in scenario 4; (b) in scenario 6. 
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